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We performed spin transport measurements on boron nitride based single layer graphene devices 
with mobilities up to 40 000 cm 2 V _1 s _1 . We could observe spin transport over lengths up to 20 fim 
at room temperature, the largest distance measured so far for graphene. Due to enhanced charge 
carrier diffusion, spin relaxation lengths are measured up to 4.5 /im. The relaxation times are similar 
to values for lower quality Si02 based devices, around 200 ps. We find that the relaxation rate is 
determined in almost equal measures by the Elliott- Yafet and D'Yakonov-Perel mechanisms. 



The potential of graphene [1] as an emerging mate- 
rial for spintronics has been established, revealing spin 
relaxation lengths A of 2 /xm at room temperature [2]. 
Spins relax over a length A = V ' D s t Si where D s is the 
spin diffusion constant and t s the spin relaxation time. 
One straightforward way to achieve spin transport over 
larger distances is to enhance D s by fabricating high mo- 
bility devices. On the other hand, r s is theoretically 
predicted to range up to hundreds of nanoseconds [3]. 
However, observations made in the recent years by ex- 
perimentalists [2, 4-13] do not match up to the high 
expectations set by theory. Measurements typically in- 
dicate t s to be in the hundred picoseconds range and 
the discrepancy between theory and experiment and the 
exact relaxation mechanism remain yet unclear. Some 
works suggest that spin relaxation is dominated by the 
Elliott- Yafet (EY) mechanism [4, 11, 14], where r s is pro- 
portional to the momentum relaxation time t p and spins 
lose their information during scattering events. Other ef- 
forts indicate that the D'Yakonov-Perel (DP) mechanism 
is stronger [3, 15, 16], where r s is inversely proportional 
to t p and spins dephase in between scattering events. 

In identifying the limiting factors on spin transport 
in graphene, the substrate deserves special attention. 
For charge transport it has already been shown that 
the standard silicon oxide substrate reduces the mobil- 
ity of charge carriers considerably [17]. The Si02 sub- 
strate is expected to also affect the spin relaxation in 
graphene through its roughness, trapped charges and 
surface phonons [16]. One approach to reduce the sub- 
strate roughness and screen impurities is to use epitaxial 
graphene on silicon carbide [12, 18]. However, the pres- 
ence of localized states is believed to affect spin trans- 
port in this system [13]. Alternatively, eliminating the 
influence of the substrate by suspending the graphene 
flake yields a 3 orders of magnitude increase in mobil- 
ity [19, 20]. Suspended spintronic graphene devices have 
been studied and a lower bound for t s of ~200 ps was 
found [10]. Determination of the actual value was how- 
ever not possible since the presence of local supports for 
the suspended device was found to dominate the extrac- 
tion of T s - 



Atomically flat hexagonal boron nitride (h-BN) was 
found to be a much better substrate than Si02 for high 
quality graphene electronics [21-23], yielding a 2 orders 
of magnitude increase in mobility. In this manuscript 
we present spin transport measurements on h-BN based 
graphene devices, which give access to a higher mobility 
regime than explored so far in graphene spin transport, 
while overcoming Si02 related issues such as the rough- 
ness and the presence of trapped charges. 

Devices arc made by transfer of a graphene flake 
(HOPG grade ZYA) onto a h-BN crystal, typically 
^20 nm thick, following a transfer recipe described in de- 
tail elsewhere [23]. The h-BN crystals are mechanically 
cleaved from a commercially available boron nitride pow- 
der (Momentive PolarTherm, Grade PT110). Electrical 
contacts to the graphene flake are made using standard 
EBL techniques. We first deposit either aluminum or 
titanium, in both cases in two steps of 0.4 nm, with an 
oxidation step after each deposition. Secondly we deposit 
65 nm cobalt in order to have spin sensitive contacts. The 
oxide barrier at the contact interface is required to tackle 
the conductivity mismatch problem, i.e. to prevent spin 
relaxation through the contacts. Typically the contact 
resistance is in the order of ~10 kS7. The standard recipe 
for high quality graphene devices on h-BN requires a fi- 
nal anneal step in Ar/H2 flow at 330 °C for 8 hours. We 
found that this treatment degrades our Co contacts so 
that they lose their spin injection and detection proper- 
ties. Therefore we omitted this step here. Also no other 
annealing steps have been used, keeping the fabrication 
process the same as for Si02 based devices. Measure- 
ments are all done in vacuum (~ 1 x 10^ 7 mbar), using 
standard AC lock-in techniques with currents up to 5 /xA. 

To characterize our devices we determine the field ef- 
fect mobility fj, of the charge carriers: \i = l/neR sq . Here 
e is the electron charge and n is the charge carrier density. 
The latter is calculated using n = —(V — Vd), where V 
is the applied gate voltage, Vu is gate voltage at which 
the charge neutrality point is found and C g ~ 67 ^F m -2 
is the geometric gate capacitance for 500 nm Si02 and 
20 nm h-BN. For device 2 shown in Fig. la we deposited 
the graphene flake partly on top of h-BN and partly on 
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FIG. 1. (Color online) (a) Optical micrograph of device 2, 
partly on SiCh (right side) and partly on h-BN (left side), 
(b) Square resistance versus charge carrier density measured 
on both a SiCh and a h-BN based parts of the same graphene 
flake. The inset shows the respective field effect mobilities. 



Si02, allowing for a direct comparison between charge 
carrier transport for both cases. The obtained values for 
R sq are presented in Fig. lb, the inset shows the calcu- 
lated mobility. The quality improvement due to the h-BN 
is directly reflected by the enhanced mobility. The inflec- 
tion point mobilities at room temperature for the three 
devices presented here are for 1: 40 000 cm 2 V~ 1 s~ 1 , 
2: 21 000 crr^V^s- 1 and 3: 14 000 cm 2 ]/"^" 1 . The 
reduction in full width at half maximum of the Dirac 
peak indicates smaller fluctuations in charge carrier den- 
sity for the h-BN supported graphene flake [19] . We find 
that despite omitting cleaning steps, the charge trans- 
port quality is well above that of SiC>2 based graphene 
spin transport devices. 

For spin transport measurements we employ the 4 
terminal non-local technique [24], schematically shown 
in Fig. 2a. We inject a spin polarized current to the 
graphene flake by sending an electrical current through 
the pair of contacts on the right side. The injected spins 
diffuse through the graphene and arrive at the detection 
circuit on the left. An external magnetic field applied 
parallel to the contacts allows for control over their mag- 
netization. Variation in contact width between 130 and 
800 nm ensures different coercivity. Spin valve measure- 
ments are taken by sweeping the parallel magnetic field 
while recording the non-local resistance R n i, as shown 
in Fig. 2b. The switches that occur in one sweep can be 
traced back to switching the magnetization of a specific 
contact. What makes this room temperature spin valve 
measurement particularly interesting is the large contact 
spacing. The total length covered is 18 /im, with mul- 
tiple switches showing up. Note that the largest length 



over which a spin signal has been observed is 20 /im for 
this device, which is the largest contact spacing that was 
available. This is a direct indication of an improved spin 
relaxation length. Note that the presence of several other 
electrodes (dashed lines in Fig. 2a) between the injection 
and detection circuits does not introduce considerable 
additional spin scattering. 

We measure Hanle spin precession by applying an 
external magnetic field perpendicular to the graphene 
plane. The diffusing spins will precess at the Larmor fre- 
quency ujl = gfisB/h, where g is the g-factor, hb is the 
Bohr magneton and h is the reduced Planck constant. 
Fig. 2c shows one set of precession data. Since in this 
case contact A contributes strongly to the spin signal, 
we measure precession for each magnetization geometry 
(at levels I to III), shown in the inset of Fig. 2c. Taking all 
three precession curves into account we can eliminate the 
contribution of the outer contact. The precession is fit 
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FIG. 2. (Color online) (a) Schematic showing the non-local 
4 terminal geometry for device 3. Contacts not used in the 
measurement are represented by dashed lines, (b) Spin valve 
measurement showing the non-local resistance versus mag- 
netic field, parallel to the electrodes, at n = 1.5 x 10 12 cm~ 2 . 
Switching of electrodes A to C shows up in the measurement, 
(c) Hanle spin precession measurement (R n i versus perpendic- 
ular magnetic field) averaged between measurements at levels 
I to III (original data in inset). The solid line shows the fit 
from which r s and D s are extracted. 
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FIG. 3. (Color online) Room temperature data extracted 
from precession measurements for device 1, 2 and 3 as func- 
tion of charge carrier density, with a respective injector de- 
tector spacing of 3.5, 3.5 and 2 /im. The barriers are made 
of AI2O3 for device 1 and 2 and of Ti02 for device 3. (a) 
Spin- and charge diffusion constants, symbols and lines re- 
spectively. The latter are extracted from charge transport, 
(b) Spin relaxation times, (c) Spin relaxation lengths. 

using the one-dimensional Bloch equation in the steady 
state regime: D s \7 2 fl s — ^ + jB x fl s = 0. Here fl s is 
the spin accumulation, D s is the spin diffusion constant 
and 7 is the gyromagnetic ratio. From a fit we acquire 
D s and r s and hence we can calculate the spin relaxation 
length A = ^/D s t s . For this particular dataset we find 
an increased D s = 0.052 m 2 /s compared to SiC>2, with 
values ~ 0.02 m 2 /s. Interestingly, t s = 390 ps is not 
much different from a SiC>2 based device. The fact that 
it is higher than the typically observed 200 ps is due to 
the use of TiC>2 barriers instead of AI2O3, resulting in a 
reduced barrier roughness [7]. 

For devices 1 to 3 the parameters obtained by fitting 
room temperature precession measurements are shown 
in Fig. 3a to c. Values for r s arc found in the range from 
50 to 480 ps for various charge carrier densities, which 
is similar to what is found for lower mobility devices on 
SiC>2. For devices 2 and 3 D s does not match well with 
the charge diffusion constant D c (solid lines in Fig. 3a), 
which is obtained using the Einstein relation and charge 
transport measurements [4] . Because our spin relaxation 
length is much larger than the contact spacing the deter- 
mination of D s is less accurate and therefore we use D c 
to calculate A. This way we obtain relaxation lengths up 
to 4.6 /Ltm. We also measured spin transport at tempera- 
tures down to 4.2 K, which only led to a minor increase 



FIG. 4. (Color online) Room temperature data obtained from 
devices 1, 2 and 3. Linear relations allow for extraction of the 
spin-orbit coupling assuming (a) the DP mechanism or (b) the 
EY mechanism. The solid lines reflect the theoretical relation 
for several spin-orbit couplings. For both mechanisms the 
data deviates from theoretical expectations, (c) The combi- 
nation of both the DP and EY mechanism allows for extrac- 
tion of the respective spin-orbit couplings using Eq.l. The 
solid lines are linear fits. 



in D s and t s . The behavior we measure for the p-doped 
part of device 3 deviates in the sense that t s decreases 
with increasing hole density. The cause for this is un- 
clear. 

In order to investigate the underlying spin relaxation 
mechanism in relation to the substrate and device quality, 
we can analyze the data from Fig. 3 in the light of the DP 
and EY mechanisms by looking at the relation between 
t s and t p . We extract the latter from charge transport 
measurements, using t p = where Vf — 10 6 m/s is 
the Fermi velocity. For DP the relation between the spin 
relaxation time t s and momentum relaxation time t p is 

given by - X Dp = (^ 4A h i F ^ t p , with Adp as the effective 
spin-orbit coupling [15]. In Fig. 4a this relation is plotted 
for three different values of Adp, together with the ex- 
perimentally obtained data. We observe that the linear 
trends described by theory are not reflected by our data. 

On the other hand we can consider an EY mechanism, in 

2 

which case the relation is given by t s ey — iff 2 , where 
Ef is the Fermi energy and Aey is the spin-orbit cou- 
pling [3, 14]. Fig. 4b shows this relation for three values of 
Aey- In this case the experimental data does show lin- 
ear trends, but none of the sets intersect with the origin. 
This cannot be attributed to broadening of the density 
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of states or finite conductivity [9] . 

Alternatively we can consider both mechanisms si- 
multaneously, with the overall scattering rate given by 
— = — 1 — , which leads to the following relation: 



- A 2 



h 2 



(1) 



When plotting '^- E - versus epr 2 for our data in Fig. 4c 
we clearly observe a linear behavior for all devices. Us- 
ing a linear relation we determine the slope and intersect 
with the y-axis, which directly gives measures for both 
Aey and Ajjp- The deviating behavior for the p-doped 
side of device 3 originates from the unexpected relation 
between t s and n (Fig. 3b) and Aey could not be ex- 
tracted. Looking at the respective relaxation rates, we 
find that tuning the charge carrier density in fact leads to 
very similar rates for both the DP and EY contribution. 
For DP we find rates in the range from 1 x 10 9 s _1 to 
2 x 10 10 s" 1 and for EY from 3 x 10 8 s^ 1 to 4 x 10 10 s -1 . 

For comparison the analysis using Eq.l is applied to 
the room temperature spin transport data achieved for 
a SiC>2 based device in Ref. [4]. In this case we find 
Adp ~ 90 [ieV and Aey ~ 2.3 mcV, which is very sim- 
ilar to the result for the h-BN based devices. Likewise, 
the relaxation rates for both mechanisms are compara- 
ble as well, with values in the order of 10 9 s _1 for both 
mechanisms. 

We measured spin transport in supported graphene de- 
vices at higher mobilities then achieved so far. Interest- 
ingly, we do not see an effect on r s , which is consistent 
with results that were recently achieved for SiC>2 based 
devices with tunable mobility [8] . The use of a h-BN sub- 
strate allowed us to exclude the influence of the SiC>2 sub- 
strate and graphene roughness on spin relaxation. These 
two types of graphene devices only have the contacts and 
the contaminants in common which therefore appear to 
dominate the spin relaxation properties of graphene de- 
vices. The contacts are known to have an effect [7], 
but the mismatch between experimental observations and 
theoretical predictions is not accounted for. We observe 
an enhanced spin relaxation time when using TiC>2 barri- 
ers, which are more smooth than AI2O3. The resistance 
of our contacts is expected to be sufficient for accurate 
measurement of non-local spin signals [5] , so the contacts 
are excluded as the dominant limiting factor for spin re- 
laxation. Other factors that can be held responsible are 
scattering by covalently bonded adsorbates and charged 
impurities [25]. An interesting approach for future ex- 
periments would be to realize clean graphene flakes by 
removing resist residues and other contaminants. First 
steps have been taken by the realization of suspended 
graphene flakes that can be current annealed [10]. For 
h-BN based devices alternative cleaning methods, as for 
example mechanical cleaning [26] , may be adopted to pre- 
serve the contacts. 



In conclusion we have fabricated graphene spin trans- 
port devices based on h-BN crystals with mobilities up 
to 40 000 cm 2 V -1 s -1 that show superior spin relaxation 
lengths at room temperature. This is directly demon- 
strated by spin- valve measurements over distances up to 
20 fjxa. and confirmed by Hanle spin precession measure- 
ments, which reveal spin relaxation lengths up to 4.5 /Ltm. 
The increase with respect to lower mobility SiC>2 based 
devices is due to an increase in the diffusion constant, the 
spin relaxation time remains virtually unchanged. This is 
an experimental indication that the substrate and rough- 
ness are not the limiting factors for spin relaxation in 
graphene, which is an important observation since most 
research has been done on SiC>2 based devices. Other 
factors may be dominant, such as contaminants on the 
graphene flake. Concerning the relaxation mechanism, 
we find that our data is best described by a combination 
of the EY and DP mechanisms. The respective relaxation 
rates are found to be very similar for both mechanisms, 
indicating that the spin relaxation is not dominated by 
a single mechanism. 
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